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PACS 72.70.+m - Noise processes and phenomena 

Abstract -A beam splitter consisting of two normal leads coupled to one end of a topological 
superconducting nanowire via double quantum dot is investigated. In this geometry, the linear 
current cross-correlations at zero temperature change signs versus the overlap between the two 
Majorana bound states hosted by the nanowire. Under symmetric bias voltages the net current 
flowing through the nanowire is noiseless. These two features highlight the fermionic nature of 
such exotic Majorana excitations though they are based on the superconductivity. Moreover, 
there exists a unique local particle-hole symmetry inherited from the self-Hermitian property of 
Majorana bound states, which is apparently scarce in other systems. We show that such particular 
symmetry can be revealed through measuring the currents under complementary bias voltages. 
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0 Introduction. — Over the last two decades, the hy- 
i_brid multiterminal structures consisting of a BCS super¬ 
conductor and two normal metal leads keep receiving ex- 
^ tensive interest from both theoretical Bi and experi- 


_mental [6||^ communities. The generic physics contained 

in this versatile platform is the interplay between co- 
l/~jherence effect in the normal leads and intrinsic coher- 
^^ence of the superconducting condensate. One of the 
^~;most appealing advantages of these structures is to act 
Os as a Cooper pair beam splitter [^[^, which splits con- 


crossed Andreev reflection (CAR), are nonlocal processes 
involving both two separated leads. At lowest order in the 
tunneling amplitudes between the normal leads and the 
superconductor, the two nonlocal processes are decoupled, 
leading to positive (negative) current cross-correlation for 
CAR (EC) process [^. A simple interpretation of this 
fact is that CAR implies instantaneous currents of the 
same sign in both leads, while EC implies instantaneous 
currents of opposite signs. This positive current cross¬ 
correlation roughly signals the Bosonic nature of Cooper 


in 


stituent spin-entangled electrons from the superconduc- pairs to some extent 16 
tor into the separate normal leads. This enables them as 


• entanglement generating sources in quantum-information 


^ processing 11 12 . It also allows for the study of elec- 
tronic Einstein-Podolsky-Rosen experiment based on cur- 
Vh rent cross-correlations 13,14. In general, the subgap 
^transport occurring in these hybrid structures with two 


normal leads includes following elementary processes [3|4| : 
an electron emitted from one of the leads is reflected back 
as a hole, or is transmitted as an electron or a hole into 
the other lead. The former one is the traditional local 
Andreev reflection (AR) occurred in one of the contacts 
between the normal leads and the superconductor, while 
the latter two, termed as elastic cotunneling (EC) and 


Recently, the interplay of AR, EC and CAR processes 
are predicted to lead remarkable transport behaviors when 
the BCS superconductor in these hybrid structures is re¬ 
placed by a InAs or InSb topological superconducting 
nanowire (TSNW) [T7[|^ , which is expected to support a 
pair of exotic end Majorana bound states (MBSs) [21-26 
that are promising candidate to realize nonlocal qubits 


free of local decoherence 27 . In particular, Zocher et al. 


studied the situation where the nanowire is sandwiched be¬ 
tween two normal leads with inserted quantum dots (QDs) 
28 . They found that the linear current cross-correlation 


at zero temperature is proportional to the overlap between 
the two end MBSs and thus vanishes for perfectly decou- 
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Fig. 1: (Color online) (a) Schematic view of our device. Two 
QDs are attached to separate normal leads and coupled to the 
nearest MBS hosted by a grounded TSNW. (b) Three elemen¬ 
tary transport processes: an electron emitted from one normal 
lead is reflected back as a hole (AR), transmitted as an electron 
(EC) or a hole (CAR) into the other lead. 


pled MBSs pair, and more particularly, is a positive def¬ 
inite quantity. This feature was also reported by Nilsson 
et al. in a similar device but without QDs 29 . However, 


such positive current cross-correlation looks to be a conse¬ 
quence of boson-like excitation, which is inconsistent with 
the fact that MBSs are fermionic excitations that obey 
canonical anti-commutation relation. 

In this paper, we investigate a beam splitter consisting 
of two normal leads coupled to one of the MBSs hosted at 
the ends of a TSNW through double noninteracting QDs, 
as shown in fig. 0 a). The inserted QDs provide a con¬ 
trollable way to achieve deep insight into the interplay of 
the elementary transport processes in virtue of the high 
tunability of QDs 30 31 . In this setup we find that. 


under appropriate bias voltages and dot levels, the cur¬ 
rent cross-correlations expected to be positive (negative) 
in conventional BCS Cooper pair splitters are indeed neg¬ 
ative (positive) in our device for weakly overlapped MBSs, 
but change signs towards strongly overlapped MBSs. Ad¬ 
ditionally, under the symmetric bias voltages the net cur¬ 
rent flowing through the nanowire is found to be noise¬ 
less. These two features on current correlations highlight 
the fermionic nature of such exotic MBSs though they are 
based on the superconductivity. Moreover, there exists a 
unique local particle-hole (p-h) symmetry in this system 
due to the self-Hermitian property of MBSs. This local 
p-h symmetry guarantees the transmission coefficient of 
charges emitting from one lead to the other lead by the EC 
process is identical to that by the CAR process through 
the reversed dot level, which can be revealed by measuring 
the currents under complementary bias voltages. 


Model. — We consider two noninteracting single-level 
QDs connected to conducting electron reservoirs and to 
one end of a grounded TSNW, as shown schematically 
in fig. 0 a). Each QD is attached to a separate normal 
lead with controllable chemical potential /ia = —eVa {a = 


L, R). Generally, a strong magnetic field is needed to drive 


the nanowire into the topological nontrivial phase 22 ■ 26 


We assume that the magnetic field is applied on the whole 
setup, and the bias voltage window is smaller than the 
large Zeeman splitting. Therefore, our considerations are 
limited to the case of spinless electrons. 

The system under consideration can be described by the 
extended Anderson Hamiltonian of the general form 


H — H\eads “b FIcentral “b FItunnel- (1) 

The first term, Flieads, describes the two noninteract¬ 
ing normal leads, iJieads = with Ha = 

J^k^kacla^^ka being the Hamiltonian of lead a. Here, 
claicka) is the creation (annihilation) operator of an elec¬ 
tron with the wave vector k in the lead a, whereas Ska 
denotes the corresponding single-particle energy. The sec¬ 
ond term of Hamiltonian (0 describes the central region 
of the system, 

.^central — ad'lda+iSM7ll2+'^Kid'l-da)ll, (2) 

a Oi 

where £„ is the discrete level of the ath QD. Since we are 
interested in the low-energy physics, we assume the in¬ 
duced superconducting gap in the nanowire is the largest 
energy scale and the applied bias voltage is controlled 
within the gap, thus the topological phase in the nanowire 
can be described by a low-energy effective model in which 
the two zero-energy MBSs hosted at two ends of the 
nanowire are represented by the Majorana operators ji 

{i = 1 , 2 ) ll 

6 ij and the self-Hermitian property 7 ^ = 'yj. Moreover, for 
a finite-length nanowire, the two MBSs would couple to 
each other by a nonzero overlap energy Em ~ 6 xp(—Z/f) 
with I being the length of the nanowire and ^ the supercon¬ 
ducting coherence length. The last term in Hamiltonian 
(0 corresponds to the tunnel coupling between the QDs 
and the nearest MBS 71 . The tunneling between the QDs 
and the normal leads reads 

l^tunnel — ^(k"a4cfea+H.C.), (3) 

koi 

with Va being the tunneling matrix elements. An elec¬ 
tron and/or hole transfer between the QD and corre¬ 
sponding lead is described by an effective tunneling rate 
Fq, = ■np\Va\^, where p is the density of states of the lead 
electrons. 

It is seen that there exists a unique local p-h symmetry 
in this system, that is, after explicit local p-h transfor¬ 
mation on operators concerning one of the QDs (say the 
right QD) as dR 4, cur 7 ^ jJ, Vr -Vr, 

£kR —SkR, and £r —>■ —£r^ when the terms involving 
the left QD keep fixed, the total Hamiltonian is invariant 
due to the nature 7 I = 7 ^ of MBSs. As we shall see, this 
local p-h symmetry can be revealed through measuring 
the currents under complementary bias voltages. 


37 , obeying the Clifford algebra { 7 z, 7 j} = 
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Transport formulae. — Following the 
nonequilibrium Green’s function formalism 38 


standard 
, the cur¬ 
rent flowing through the ath lead can be calculated from 
evolution of the total number operator of electrons in that 
lead, Na = For superconducting system, it 

is convenient to define the GFs in the Nambu space in a 


((da, 4)) ((da, d/3)) 


with 


a = (3 {a ^ 13) representing the local (nonlocal) GFs be¬ 
tween the same (different) QD. Accordingly, we denote 
the two elements in the upper row of Gap as Gap-ee and 
Gap-,eh while those in the lower row as Gap-,he and Gap-,hh- 
Within such representation, the current flowing through 
lead a reads 


Ia = 


2 h 


where a = 


0 


Sap = SapiO) = F/da;Tr{(5„^(G<„crE>cr-f crE<crG>„) 

+G4o’[E/3G^qEq,]^(t — (T[EQ,GQ/3]^cr[E,gG^Q,]^ 

— [GQ,/3E^]^cr[G/3QEQ,]^(T -|- (t[Eq,Gq,^E^]^(tG 4 }- 


of the electron and hole states in the ath normal lead. 
Note that due to the large gap assumption in our low- 
energy effective model, the quasi-particle tunneling be¬ 
tween the normal leads and the TSNW is forbidden in 
this subgap regime. Differently from the situation for cur¬ 
rent, it is, in general, impossible to combine all the terms 
in current-noise power spectrum to module squares of GFs 
, as well as of scattering matrixes 


5p^ . The current 


formula eq. § reduces to instructive form when applying 
specific bias voltage configuration. 




E 

Ti—e.h ‘ 


dw(rjr + r„"r)(/a-/a), (8) 


dwTr{CT[GaaEa]^ - Cr[E„Gaa]<}, (4) 


for the symmetric bias voltages (SBV) Vl = Vr 
where the EC process is completely blockaded, and 


= V 


la = 


accounts for the different type 


q 

0 

charge carried by the electrons ( 9 ® = e) and the holes 
{q^ = —e). The correlation function between current 
fluctuations in lead a and j3 is defined as Sapit) = 
{{dla{t),61p{0)}. Its Fourier transform gives the current- 
noise power spectrum as Sap{i^) = f Sap{t). The 
general form of the power spectrum in the zero-frequency 
limit reads 


2 h ^ 

m—e,h ‘ 


MTz^+T:r){ra-fa). ( 9 ) 


(5) 


for the antisymmetric bias voltages (ASBV) Vr = —Vr = 
V where instead the CAR process disappears. 

In the noninteracting case the exact GFs can be accessed 
with analytical expressions, and so do the transmission co¬ 
efficients. The matrix elements of local and nonlocal GFs 
can be obtained through the Dyson equation Gap-,mn = 
3ap3mn9cia-,mra “t" 9cta-,7nmVa'-^-,mlG^^-,\\V^p-^\ri9pp-,nn^ where 
we have employed two facts that (i) the bare GFs for 
QD a decoupled from the MBSs is in a diagonal form 

0 

1 


as 5aa(w) = 


o; —eaiirQ 

0 


The enclosed GFs in eq. 0 and eq. 0 can be 
expanded using the so-called analytic continuation 
rules 38 . Noting that such an 

y)mm = ^p=L,R^n=e 

the Keldysh equation G<„ 


and more im- 

portantly (ii) the particular hopping matrix between QD 

A„ 0 


a and the MBSs is Va-, = 


(GL 

SS). 


- G: 

.G 


otct/mm 

and 

GL 


useful identity 

(ypr _ 

-e,h 1*^/3 


and V-ya = 


—Act 0 

which accounts for the self-Hermitian property of 
GFs of MBSs G,y-y can be cal- 


ctcf.mm 


y< 

3 \mn^ 6 cf,nm' 


In the noninteract- 


^^P—L,R ^^n—e,h ' 

ing case, the self-energies results from tunnel coupling 
between the QDs and the normal leads have the form 

J J ) and E< = 2zTa f 

ter some algebra we arrive at the final current formula 
consisting of three elementary transport processes, shown 
schematically in fig. 0 b) , as follow 


9'1't 


yVaa - 


= T*Fo 


Af- 


1 / 1 ' 
a-y- 

MBSs. The dressed 
culated through G;^ 
uj — E..),(w) —isM 
iSM w 

^^a —L,i3 ( Oj —iira OJ-t-Caizra 

we immediately obtain the matrix elements of GFs as 


'l2a=L,R ^ia9a 
with the self-energy E!)’“(a;) = 

-). Gombining these terms 


Gap-,mn 3ap3mn9aa-,mm 3~ 


Va^-,ml9i 


Oiy‘,mLyacx;mm 


;nn 


UJ — 


M 


1^ 


( 10 ) 


lo. = ^m=e,h £ I 

AR 

+uTifa - /D+ r:rifa - m, ( 6 ) 

EC CAR 

with corresponding transmission coefficients 

r„X(a;)=4r„F^|G^^._(cc)4 (7) 

Here d and m denote the index opposite to a and m. 
fa^{uj) = [exp( ‘^^y° ) -1-1]“^ are the Fermi distributions 


Results and discussions. — We investigate the cur¬ 
rent and current correlations in such topological supercon¬ 
ducting beam splitter under complementary bias voltages. 
In what follows, we consider the case F^ = Ffl/ = F and 
take F as the energy unit. Furthermore, to be more realis¬ 
tic we adopt asymmetric couplings between the two QDs 
and the MBS 71 in calculations. All results are obtained 
in the zero temperature limit. 

In fig.we show the current Ir flowing through the left 
lead as functions of the dot levels under complementary 
bias voltages. In fig.j^a) and fig.j^b), the currents with 


P-3 


















Zhan Caol Tie-Feng Fangl Qiao Chen2 Hong-Gang Luol,3 



Fig. 2: (Color online) Density plots of the current II (in unit 
of eF/ h) flowing through the left lead as functions of the dot 
levels under (a) symmetric bias voltage Vl = Vr = 10, (b) 
antisymmetric bias voltage Vl = —Vr = 10, (c) Vl = 10, Vr = 
0 and (d) Vl = 10 , Vr = 4. Other parameters are: Sm = 0, 
Al = 4 and Xr = 3. Note the symmetries of currents displayed 
in (a)-(c) persist for nonzero em since the inherent local p-h 
symmetry is independent of the overlap between MBSs. 


dot levels {sl,£r) are identical to those with (—£l, —£«)• 
This symmetry is attributed to the combination of the 
global p-h symmetry of the system and the particular bias 
voltage configuration. To be more precise, after the global 
p-h transformation the Hamiltonian is invariant but the 
roles of electrons and holes are exchanged. As a conse¬ 
quence, the transmission coefficients submit to the iden¬ 
tity that 7 ^’]^"(£l,£r) = 72^”(-£l,-£_r), which can also 
be readily inspected by eq. ([^ and (10 1 . This identity 
guarantees the invariance of currents addressed by eq. (|^ 
and eq. (|^ after reversing both dot levels. Another inter¬ 
esting symmetry of the currents is that the currents with 
dot levels (£L,£fl) in hg- Ha) are identical to those with 
{£l,-£r) in fig. Hb)- However, as we noted above, the 
involved elementary transport processes under SBV and 
ASBV are essentially different. This symmetry of currents 
is due to the nontrivial local p-h symmetry of the sys¬ 
tem, which indicates that T^J^{£l,£r) = T£^{£l, —£r)- 
Therefore, comparing and ®, it renders the contribu¬ 
tion to the total current from the AR and CAR processes 
under SBV with dot levels {£l,£r) being the same with 
that from the AR and EC processes under ASBV with 
dot levels {er, —£r)- In fig.^c), the chemical potential of 
right lead is set to be aline with the Fermi level. In this 
case, the currents with dot levels {£l,£r) are identical to 
those with {er, —er), which is also introduced by the local 
p-h symmetry of Hamiltonian 0 , and the particular bias 
scheme. However, as shown in fig. Hd), the symmetries of 
currents can not be uncovered under general asymmetric 
bias voltages. It is worth noticing that these intriguing 


symmetries of currents originated from the local p-h sym¬ 
metry of the system is apparently not presented in other 
non-MBSs systems. 

The maximum currents under SBV are along the line 
Er = —Er [fig.Ha)] while under ASBV are along the line 
Er = Er [fig.lWb)]. In the former case, where the EC 
process is prohibited, an electron transferred through the 
left dot level er will be resonant transmitted as a hole to 
the other lead when the right dot level equals —er since 
the Cooper pairs are residing near the Fermi level. On 
the contrary, in the latter case the EC process is allowed 
while the CAR process is instead prohibited. An electron 
injected through the left QD is to be resonant transmitted 
as an electron when er = er. Besides, the currents in each 
plots in fig.also exhibit resonant peaks centered around 
Er = 0, which are results of the local AR process occurring 
between the left lead and the TSNW. 



Fig. 3: (Color online) (a) The II as a function of bias voltage 
under SBV (red) and ASBV (blue). The solid and dashed lines 
denote em = 0 and em = 2, respectively, (b) The correspond¬ 
ing differential conductances dlR/dV. Inset: Evolution of the 
linear response conductance under SBV with increased over¬ 
lap as Em = 0,0.03,0.1,0.2,0.4. The parameters are el = 4, 
Er = —2, Al = 4 and Xr = 2. 


The I-V characteristics under SBV and ASBV with and 
without overlap between the two end MBSs are displayed 
in fig.[^ In the low-bias regime, the currents increase lin¬ 
early as the voltage, as shown in fig.^a). The slopes corre¬ 
sponding to strong overlap Em ~ Al, Alj are much smaller 
than those to zero overlap, which are clearly shown by 
the differential conductances in fig.[^b). For zero over¬ 
lap, there are distinct zero-bias conductance peaks arising 
from the resonant EC, AR and CAR processes via the 
zero-energy MBS 71 with different weights depending on 
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the specific bias voltages and dot levels. However, as the 
overlap becomes nonzero, the resonant peak splits into 
double peaks around the Fermi level symmetrically with a 
central dip [inset of fig. ib)]. In the nonlinear regime, all 
the differential conductances peaks under SB (red lines) 
are higher than those under ASB (blue lines), which is 
due to the dot levels we adopt have opposite signs that 
favors the CAR process. If the right dot level is reversed, 
the I-V curves under SBV and ASBV will exchange to 
each other, as we discussed in fig.[^ Around the zero- 
bias, the differential conductances under SBV and ASBV 
are nearly identical [hg.[^b)]. This feature is ascribed to 
the zero-energy transmission coefficient of EC process is 
identical to that of CAR process as 


r„’S"*(o) = 7;™’^(o) = ^fs:FWTWW). (n) 


■^bh \/ 72 S™( 0 ) — Y-L>^i(e%+T\)+TRXj^{el+Tl) “ b). 

Note this equivalence between zero-energy transmission 
coefficients does not require the sign reversal of either dot 
level. We would like to mention two nontrivial features on 


equation (11). (i) If A/, = A^, the zero-energy transmis¬ 


sion coefficients are independent of the QDs-MBS coupling 
strengthes since Xl and Xr cancel out each other, (ii) At 
the high symmetric points A± (F^ = F/j = F, = ±e_R, 
Al = Xr), all three transmission coefficients turn out to 
be 1/4, which is independent of the particular dot levels 
as well as the QDs-MBS coupling strengthes. 

In hg.|^ the current cross-correlations under SBV and 
ASBV are calculated. Similar to the currents in fig.[^ the 
cross-correlations under SBV with dot levels er = —sr are 
identical to those under ASBV with dot levels £l = £r ex¬ 
cept a negative sign, which also resultes from the local p-h 
symmetry of this system. As the voltage increases, more 
and more tunneling channels are opened but the corre¬ 
lations vary nonmonotonically with V, depending on the 
cooperative or competitive relationship between these tun¬ 
neling channels through which charges transferred. The 
cross-correlation between same type charges is always neg¬ 
ative while between opposite type charges is always pos¬ 
itive 41 . As a result, the total cross-correlation can be 


either positive or negative up to the relative strengths of 
the above two terms. In the linear regime, where only the 
MBSs related low-energy channel is activated, the cross¬ 
correlation under SBV is negative for small overlap £m 
(see red solid and dashed lines). This is a counterintuitive 
feature since the dot levels are chosen to be er = —er 
to evoke resonant CAR process that likely to contribute 
positive cross-correlation [^. Nevertheless, as the over¬ 
lap Em is increased the linear cross-correlation becomes 
positive as one may expected in the conventional Cooper 
pair beam splitters (see red dotted line). Similarly, the 
linear cross-correlation under ASBV with resonance con¬ 
dition Er = Er for EC process is unexpectedly positive for 
weak overlap but turns to be positive for strong overlap 
(see blue lines). These seemingly anomalous signs of linear 
cross-correlations are appreciable since the MBSs are pre¬ 



Fig. 4: (Color online) Current cross-correlations Srr under 
SBV with El = —er = 2 (red) and ASBV with er = er — 
— 2 (blue) as a function of the voltage with different overlaps 
between the two MBSs as £m = 0 (solid lines), em = 2 (dashed 
lines) and £m = 4 (dotted lines). Inset: Differential current 
correlations defined in Eq. (12 1 versus overlap strength at V = 
0.01. The QDs-MBSs couplings are An = 4 and Xr = 2. 


dieted to be fermionic excitations, however, the strongly 
overlapped MBSs are akin to standard Andreev bound 
states. We note these sign reversals of cross-correlations 
versus the overlap are absent when the TSNW is sand¬ 
wiched between two normal leads with or without inserted 
QDs 


28 29 


To gain more analytical insights on the linear current 
correlations, we obtain in the zero-bias limit and Em = 0 
the differential correlations Wq /3 = dSap / dV as follows (in 
unit of e^/h) 

Wrr{V ^ 0+) = -Wrr{V ^ 0+) = P, (12) 
for the SBV, and 

WllCC^ 0+) = WLfl(V^0+) =P, (13) 

for the ASBV. Interestingly, P = 4:TETi'^) = 472^™ (0), 
which can be tested by tuning the dot levels according to 


eq. (11) to extremely suppress the AR process and com¬ 


paring the conductance and differential current correla¬ 
tions. Furthermore, the linear cross-correlation is always 
negative (positive) under SBV (ASBV), irrespective of the 


specihe model parameters. For SBV, eq. (12) indicates the 


linear self-correlation Srr = WrrAV equals the cross¬ 
correlation Srr = WrrAV except a minus sign, which is 
robust at small overlap [inset of fig.|^ Considering the cur¬ 
rent conservation condition Ir + Ir + Iw = 0, this feature 
implies the net current flowing through the nanowire is 
noiseless since S^w = —Srr — Srr — 2Srr — 0, which un¬ 
derlies the fermionic nature [4 42 of MBSs though they 
are based on the superconductivity. Such noiseless fea¬ 
ture goes beyond the previous noiseless current due to 
the MBS-induced resonant Andreev reflection 


17 


since 


Eq. (12) is always valid regardless of the on-resonance 
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(T^” = 1) or off-resonance ^ 1) situation. Par¬ 

ticularly, this feature is distinct from that the difference 
of currents II — Ir is noiseless in a BCS Cooper pair beam 
splitter with highly transmitting interfaces . At the high 
symmetric points A±, P reduces to an universal value 1 
which is similar to the one reported by Liu et al. who 


justified it as a unique signature of MBSs 43 


Summary. — In conclusion, the currents and current 
correlations in a normal lead-double QD-TSNW hybrid 
beam splitter are investigated. We find the linear cross¬ 
correlations at zero temperature change signs versus the 
overlap between the two MBSs. Under symmetric bias 
voltages, the net current flowing through the nanowire is 
noiseless. These two features highlight the fermionic na¬ 
ture of such exotic Majorana excitations. Moreover, we 
address a unique local p-h symmetry inherited from the 
self-Hermitian property of MBSs, which can be revealed 
through measuring the currents under complementary bias 
voltages. All these predictions are observable in transport 
measurements by applying appropriate bias voltages and 
carefully tuning the discrete levels in QDs. Additionally, 
we note the nonlocal advantage of such beam pair splitter 
device facilitates experimentalists to exclude possible local 
mist in the unambiguous identification of MBSs. 


* * * 
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